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Abstract 

We consider the electronic transport in multi-terminal mesoscopic networks of weakly disordered metallic 
wires. After a brief description of the classical transport, we analyze the weak localisation (WL) correction 
to the four-terminal resistances, which involves an integration of the Cooperon over the wires with proper 
weights. We provide an interpretation of these weights in terms of classical transport properties. We 
illustrate the formalism on examples and show that weak localisation to four-terminal conductances may 
become large in some situations. In a second part, we study the correlations of four-terminal resistances and 
show that integration of Diffuson and Cooperon inside the network involves the same weights as the WL. 
The formulae are applied to multiconnected wire geometries. 
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1. Introduction 

Classical laws of transport on electrical networks 
have been established by Custav Kirchhoff in 1845. 
They rely on three fundamental hypotheses. Two of 
them are energy and charge conservation leading re¬ 
spectively to the so-called voltage and current laws. 
The third one is Ohm’s law which states that the 
current along a wire [pv) of the network is pro¬ 
portional to the voltage drop between the vertices 
p and V connected by this wire : 

V = (1) 

being the length of the wire, s its cross section 
and (To the Drude conductivity. We now under¬ 
stand that this third hypothesis relies on the as¬ 
sumption of diffusive motion of the charge carri¬ 
ers at the microscopic level, within a classical de¬ 
scription. In particular, it is not appropriate to de¬ 
scribe quantum effects (like Aharanov-Bohm effect) 
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or non-diffusive regime (ballistic regime, quantum 
Hall effect). After the fundamental breakthrough 
proposed by Rolf Landauer to describe the elec¬ 
trical conductance as a transmission coefhcient, a 
generalization of laws of transport beyond classical 
transport was highly desirable. It has been formal¬ 
ized in a beautiful work by Markus Biittiker (for 
reviews, see Refs. [211 [2711^149]). 

If we restrict ourselves to the regime of linear 
transport, a convenient description is to start by 
introducing the conductance matrix G which relates 
the currents at the contacts (also called terminals) 
of the circuit to the values of the voltage at these 
terminals 

/a = X ■ (2) 

P 

This relation is completely general in the limit of 
linear transport. The two first Kirchhoff’s laws im¬ 
ply : 

Y.G^0 = Y,Gc.s = O, (3) 

a p 

expressing the conservation of current, and the in¬ 
variance of the current distribution against a global 
shift of the potentials (gauge invariance). As a gen- 
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eralization of Landauer’s formula, elements of the 
conductance matrix are related to transmission co¬ 
efficients 

Ga /3 = for /3 ^ a , (4) 

where the factor 2s stands for the spin degeneracy. 
At this level, the formalism is completely general 
and no hypothesis is made on the nature of electri¬ 
cal transport which is totally encoded in the trans¬ 
mission coefficients These coefficients are re¬ 
lated to the scattering matrix which may be deter¬ 
mined explicitly within specific models. For exam¬ 
ple, in the regime of the integer quantum Hall effect, 
the current is carried by edge states which makes 
the problem effectively one-dimensional and allows 
for a simple construction of the scattering matrix 
[iziiin]- More generally, the scattering matrix may 
be constructed efficiently by assuming strictly one¬ 
dimensional character, like in a network of strictly 
one-dimensional (ID) wires [59]. Another powerful 
approach applies to devices in which the electron 
dynamics can be considered as ergodic, leading to 
a random matrix formulation of quantum scatter¬ 
ing [12]. In the present article, we consider the 
case of metallic samples made of weakly disordered 
wires such that the electron dynamics is diffusive^ 
as it is the case in narrow metallic wires deposited 
on a substrate [72] (see also the recent experiments 
[53] 1. or wires etched in a two-dimensional-electron 
gas [3211331 Eg. 



deal with voltage measurements : current is in¬ 
jected and collected at two specific contacts and 
voltages are measured at any pair on contacts play¬ 
ing the role of voltage probes (Fig. [^. Therefore, 
the relevant quantities characterizing the response 
of the device are the four-terminal resistances de¬ 
fined as 


= 7 

V^-V, I “ 

= - J - With -I 

[/;, = 0 VA a,/3 

(5) 

By appropriate inversion of relation ([^, Biittiker 
could relate the four-terminal resistances to the el¬ 
ements of the conductance matrix, therefore to the 
transmission coefficients mils] : 


TZn 




^ /x/3-^ I/Q: 


2.e2 


V 


( 6 ) 


where V is any minor of the dimensionless conduc¬ 
tance matrix. The expression is valid when all in¬ 
dices are different. 

In this paper, we are interested in transport prop¬ 
erties of mesoscopic diffusive wires, where quantum 
interferences lead to small deviations to Ohm’s law. 
This is the so-called weak-localisation regime. For 
classical transport, simple application of Ohm’s law 
leads to the expression of the transmission coeffi¬ 
cients, which amounts to classical addition of re¬ 
sistances and conductances (Kirchhoff). They are 
expressed in terms of the elements of a matrix which 
encodes the conductances of all the links of the net¬ 


work, defined below in Eq.(ll|. 

For a single wire of length L, it is well-known 
that the weak-localisation correction to the classical 
transmission coefficient can be written as [2] 


AT = -— / dxPc{x,x). 


(7) 


Figure 1: The four-terminal resistance is the ratio 

of the voltage between two contacts {i and u and the current 
injected at contact ol and collected at contact p. 

Biittiker emphasized the importance of the mea¬ 
surement process when a quantum circuit is con- 
nected to the outside macroscopic world [TB]. The 
measured resistance is not only a property of the 
system itself but also depends on the way it is con¬ 
nected to the outside world. Moreover, although 
the concept of conductance is natural from a theo¬ 
retical point of view, experiments most frequently 


where AT = (T) — Tdass is an average over disor¬ 
der configurations. The so-called Cooperon Pc{x, x) 
measures the contribution of interfering closed elec¬ 
tronic diffusive trajectories. We have shown that 
in a network of diffusive wires, this simple relation 
generalizes to [60] 

O _ QrpclaSS j. 

ATgp = ^ X! fi / ^xPc{x,x), (8) 

?loc . Oli Jwirep) 

where i labels all the wires of the network, ^joc = 


2 









UdNcf-e is the localisation length in the infinitely 
long wire with conducting channels [12], if. the 
elastic mean free path and ad a dimensionless pa¬ 
rameter of order unity, which will be given below. 

From the knowledge of the quantum corrections 
([^, we will show in this paper that the quantum 
correction to the classical four-terminal resistance 
is 


• A fundamental aspect of quantum transport 
concerns the symmetry of transport coeffi¬ 
cients : symmetry with respect to current 
sources and voltage probes exchange, and sym¬ 
metry with respect to the reversal of the mag¬ 
netic field B. Extending these ideas to coher¬ 
ent conductors, Biittiker has demonstrated the 
symmetry relation for the 4TRs [HdS] 






o class 

6oc ^ dli 


dxPc{x,x). (9) 


/ wire (f) 


This expression is quite simple since the weights at¬ 
tached to each wire have a simple interpretation : 
they express the sensitivity of the classical four- 
terminal resistance when the resistance of this wire 
is modified. 

Similarly we have found convenient expressions 
for the correlation functions of the transmission 
coefficients, from which one can obtain the cor¬ 
relation functions of the four-terminal conduc¬ 
tances. Like the weak-localisation correction, these 
expressions involve contribution of all the wires, 
which are weighted by similar factors k, 

Eqs. ( |48|49|50|51D . These equations, which are, 
with Eq. If , the main results of the article, will be 
illustrated by several examples in simple devices. 


Before going specifically to the analysis of quan¬ 
tum transport in networks of quasi one-dimensional 
weakly disordered wires, we close the section with 
some general remarks as the concept of four- 
terminal resistance (4TR), of which Biittiker has 
been one of the main promoters, has been extremely 
fruitful in mesoscopic physics. Let us mention few 
directions : 


• In the early developments of the Landauer- 
Biittiker approach, the concept of four- 
terminal resistance has helped clarifying the 
question of contact resistance. The concept 
of contact resistance is the “mesoscopic ver¬ 
sion” of the electric resistance appearing when 
the electronic fluid is injected from a macro¬ 
scopic conductor into a small hole (known as 
Sharvin resistance [M])- The role of contact 
resistances has been nicely explained in sev¬ 
eral papers of Biittiker [sainKiiiiig] (see also 
chapter 5 of [40]). 


^Note that the perturbative approach is valid for 
min(L,^,size) < ^loc- 


7^a/3.^.(-S)=7^^..a/3(^)• (10) 

• The concept of 4TR provides an illuminat¬ 
ing description of the integer quantum Hall 
effect dZl from the edge state picture intro¬ 
duced by Halperin |3S| , as it allows to compute 
straightforwardly the longitudinal and Hall re¬ 
sistances (see also Buttiker’s beautiful review 
article [20]). Furthermore, this framework per¬ 
mits to analyze in simple terms other more sub- 
btle effects, like the scattering between edge 
states at opposite boundaries due to a con¬ 
striction, leading to so-called anomalous Hall 
effect |44l [TO] [2^ , or by impurities [T8| , the de¬ 
scription of transport in a Hall cross [34], etc. 

• Another issue which has been put forward by 
Biittiker with others, and which will be cen¬ 
tral in the present article, concerns the non¬ 
local nature of quantum transport, and the 
influence of voltage probes on the transport 
properties of a coherent conductor. Moti¬ 
vated by a set of experiments in multicon- 
nected metallic wires [m US] (see also Whas- 
burn and Webb’s review [72]), various au¬ 
thors have analyzed the role of voltage probes 
in devices made of disordered wires by vari¬ 
ous approaches : Maekawa, Isawa and Ebi- 
sawa [3H]) Biittiker m, Divincenzo , Kane 
and Lee [32 HI]) Chandrasekhar, Prober and 
Santhanam [smsaES], Hershfield and Am- 
begaokar [3H1I3Z]- A more general discussion 
of nonlocality of weak localisation in networks 
of metallic wires was made possible within the 
theory developed by us in Ref. EH]. Our for¬ 
malism has allowed us to study how Altshuler- 
Aronov-Spivak oscillations of the magnetocon¬ 
ductance are affected by the network geometry 

[5011501 [5n 155110511551155] (this was beauti¬ 
fully demonstrated by earlier experiments in 
arrays of lithium rings by Bishop, Dolan and 
Licini m), or the role of electronic interac¬ 
tions [IIIISIIEIIIMIIOIIMIISTIIMIIIO]- 
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The outline of the paper is as follows : in the next 
section we introduce a classical description of elec¬ 
tronic transport in networks of metallic wires. Sec- 
tionj^is devoted to the analysis of the weak localisa¬ 
tion correction, as a warm up exercice preparing the 
more complicate study of four-terminal resistance 
correlations presented in Section The main for¬ 
mulae, Eqs. (|^, etc., are illustrated on simple cases. 
Section closes the paper with some concluding re¬ 
marks. 


2. Classical transport in networks 

2.1. Electrostatic potential 

A typical network is represented in Fig. In 
this section, we introduce a specific notation in or¬ 
der to distinguish between internal vertices and ver¬ 
tices corresponding to reservoirs, by labelling these 
latter with a prime. Before going to the discus- 



Figure 2: A four-terminal mesoscopic network of metallic 
wires. Vertices are marked with magenta dots. Internal 
vertices are characterised by Xfi = 0. Primed vertices corre¬ 
spond to reservoirs for which = oo. 


Fermi energy. The current in the wire in' is given 
by 0 , hence current conservation at each “inter¬ 
nal” vertex neighbour of may be 

rewritten as i-^o)= 0, where the matrix 
is defined as 



is the adjacency matrix element, equals to 1 
if a wire connects the two vertices and 0 other¬ 


wise ; thus it constraints the sum in (111 to run 


over vertices neighbours of /i. The parameters 
have been introduced for convenience for the follow¬ 
ing and describe connection to reservoirs (A^ = 0 
for an internal vertex and —>■ oo for a ver¬ 

tex in a reservoir, cf. Fig. [^. Up to a factor 
sao, the matrix Aio thus simply gathers all the 
wire conductances Oqs/1^,j, where s is the cross 
section of the wires. If we split the vector gath¬ 
ering the electrostatic potentials at the vertices 
into two parts related to internal and external ver¬ 
tices, (UinIVres) = (• • • , 14, • • • | , 14/, • • • ), we 

may write = - {Mo) i.e. [] 

K/ = Eres./3' ^^ere the sum 
runs over the reservoirs. It will be convenient for 
the following to introduce the notation Pj^{a,j3) = 
{Mo^)^p which represents the so-called Diffuson, 
measured at the two vertices : 


V^= V Pd{a,fi) — Vp, . (12) 


The Diffuson is solution of the diffusion equation 
— dlPd{x,x') = 5{x — x') (13) 


sion of quantum transport, let us analyze the clas¬ 
sical transport in the network. For this purpose we 
start by solving the Poisson equation AV{r) = 0 
for the electrostatic potential V{r) inside the net¬ 
work. Boundary conditions are V{r) = 14' for r £ 
reservoir a'. Introducing the coordinate x measur¬ 
ing the distance along the wire piv (from fi to v), the 
potential varies linearly asV{x) = (1 — ) -I- 

14 xjl^^j. The (classical) current density is given by 
Tick’s law j{f) = —eDV5n{r), where D is the dif¬ 
fusion constant. The density in excess 5n is related 
to the potential through the effective (screened) 
Coulomb interaction eV{r) = {1/vo) 5n{r), leading 
to jif) = —(7o^V{r), where ctq = e^voD is the 
Drude conductivity and vq the density of states at 


with Dirichlet boundary conditions at the reser¬ 
voirs : Pd{x,a') = Pd{a',x) = 0 for all reser¬ 
voirs a' (for details, cf. Appendix of Ref. [55]1. 
Using the linearity of the Diffuson on the wires 
and that it vanishes at the reservoirs, we write 
Pdi»,x) = {x/lpp:) Pd[»,fi) for X e P'P (the reser¬ 
voir is at a: = 0). Hence, we can rewrite the coeffi¬ 
cients in ( |l^ as Pd{ce, I3 )/Ii3P' = d^'Pd{oi,x'), where 
x' is any position in the wire fi'(3. In the following 
we will prefer to write the slope of the Diffuson on 


^since A^/ —^ oo, the inverse of the matrix is simply the 
inverse of the block (AIo)j^ related to internal vertices. 

All other matrix elements are zero, (Ad^^)^,^ = 0 for any 
reservoir p,'. 
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the wire /3/?' as Pd{a, P')/£e, where Pd{»,P') de¬ 
notes that the argument is taken at a distance £e of 
the reservoir /3', ig being the elastic mean free path, 
i.e. the smallest length scale of the problem. The 
Diffuson at a distance ig naturally appears in the 
diagrammatic calculation of the classical transport 
coefficients EO]. The relation with this reference’s 
notations will be clear by using below this expres¬ 
sion for the Diffuson’s slope. Since the vertex a 
may be any point of the network, we may rewrite 


(121 as 


res. /3' ^ 


(14) 


the expression (14) is only valid when x is at dis¬ 


tance larger than £g from the reservoirs. 

3.3. Current distribution and generalised conduc¬ 
tances 

Integrating the current density jx{x) = 
—(jQdxV{x) across the section s of the wire we 
obtain the current under the form : 




(15) 


res. 

where we have introduced 


?wire /i^',res. 


CToS dPdix,pJ 
£g dx 


for X G fin . 

(16) 

The quantity Gfiu,/ 3 ' is a generalised conductance 
matrix relating the external potentials to the inter¬ 
nal currents. They obviously satisfy the symmetry 
property = -Gvti.p'- Although the physical 

interpretation was not provided in Ref. [60] , the ex¬ 


plicit expression of (16) in terms of the matrix AIq 
was given : 

G^v,li' _ (-^0 ~ (-^0 Ippi 


CToS 


^/9/3' l^u 


(17) 

If we consider the case of a wire connected to a 


reservoir, fiv —)■ a'a, Eq. (151 coincides with the 


usual relation between currents and voltage in the 
terminals : 


res. j3' 


G class 
oc'l3' 


Vp. . 


(18) 


The conductance matrix is obviously related to the 
generalised conductances by 


The expression of the classical conductance matrix 
in terms of the matrix A4o may be deduced by set¬ 
ting fiv -G a'a in 0 : we recover the expression 

of Ref. EDI 




h 


adNgig 


(•^o')a/3 
laa' ^j3p' 


( 20 ) 


We have used aos = {2^0^/h)adNc£g , where 2^ is 
the spin degeneracy, Ng the number of conducting 
channels and Ud = VdjVd-i involves the volume of 
the d-dimensional sphere of unit radius (thus ai = 
2, Q !2 = 7’‘/2 and as = 4/3). 


Example. As a simple illustration of Eq. ( |20[ ), we 
consider the ring of Eig. From the definition 
(11), we write the internal part of the matrix AIq 
( i.e. the block related to vertices 1 and 2) : 

^ V -lAc||d l/4||d + l/4y ’ 

( 21 ) 

( 22 ) 


where l//c||d = 1/^c + 1/^d- Eq. ( [^ leads to 




1 


H” ^c\\d 

giving the expected conductance. 


it 


o 


Figure 3: A ring connected to two reservoirs. 


3.3. Four-terminal resistances 

Using the 4TRs can be deduced from the 
conductance matrix (20) (see also the discussion in 
Appendix A). For simple enough networks, the de¬ 


termination of the resistances is however more sim¬ 
ple than that of the conductance matrix and does 
not require the knowledge of the latter, as the sim¬ 
ple example analyzed just before has shown. 

3. Weak localisation 


G class _ ^ 

a'j3' — yoc'a,(3' • 


(19) 
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Weak localisation is a small quantum correction 
to transport coefficients originating from quantum 























interferences between time reversed electronic tra¬ 
jectories 0 |2]. The main interest in this small 
quantum correction to transport coefficients is that 
it gives a measure of the phase coherence length L^, 
the fundamental characterisitc length scale which 
sets the boundary between quantum and classical 
physics. It is worth stressing that there is no in¬ 
trinsic definition of L^p, which can only be obtained 
by extracting a characteristic length scale from the 
analysis of a physical quantity sensitive to quan¬ 
tum interference, such as weak localisation. A pre¬ 
cise experimental determination of thus requires 
a perfect knowledge of the functional form of the 
transport coefficients as a function of the various 
length scales, the magnetic field, etc. 


Expression (23) is valid for o' ^ /3'. The di¬ 
rect diagrammatic calculation of the reflection 
probability is more difficult and involves a de¬ 
scription of the matching between the metallic 
system and the contacts which goes beyond the 
derivation of 


(23) 


However the WL 


correction to the diagonal conductance matrix 
elements can always be deduced by using cur¬ 
rent conservation Gap = 0. 


• Interestingly, we see that the contribution of 
each wire is weighted by the “internal conduc¬ 
tances” introduced above : 


P' ~ 7 ^ / Qiiu.a'Q 

(.■SCoj 


(25) 


3.1. Conductance matrix 


The weak localisation correction to the conduc¬ 
tance matrix elements Ga'p' = —{,‘^sG^/h)Ta'p> is 
given by m 


^Ta'P' = — 






e Network 


dx 


dx 


where Pc{x,x') is the Cooperon, solution of 


(23) 


{l/Ll-Dl)P,{x,x')=6{x-x'), (24) 


Dx = dx — 2ieA(a;) being the covariant derivative 
and the phase coherence length. Few remarks : 


X 



da; Pc(x, x) . 
iuG 


This expression shows that a uniform integra¬ 
tion of the Cooperon in the network is possible 
only if the distribution of the classical currents 
in the wires is uniform. 


More conveniently, we showed in Ref. [SO] that 
these weights are related to the derivatives of the 
classical conductance matrix, Eq. ([^, or equiva¬ 
lently : 

9 ^^class p 

^Gap = ^ X! if / da;Pc(a;,x), ( 26 ) 

sloe ^ O Li J 


• The effect of the magnetic field is twofold : 
(z) in the presence of loops in the network, 
the Cooperon depends on the magnetic fluxes, 
which leads to Altshuler-Aronov-Spivak oscil¬ 
lations [H 0 HO] (see also [SllEH])- (**) The 
penetration of the magnetic field in a nar¬ 
row wire of width w can be accounted for 
through the substitution 1/T^ —t 1/T^ -I- 
1/Lg, where the magnetic length is Lg = 
[•\/3/(27r)] (j)o/{Bw) and (j)o = hje the quantum 
flux 0. 


The expression (23) is of great generality : 
it is not only valid for a system made of 
quasi-ID wires (network) but only assume 
that the contacts have a quasi-ID geome¬ 
try. In such a more general situation, the 
derivatives should simply be replaced by gra¬ 
dients dxPd{Qi,x)dxPd{x,§^ -L- VPdi^.f) ■ 
VPd(r,/I'). 


where ^loc = adN^ie is the localisation length in 
the infinitely long wire. From now, we will drop 
the prime on the vertices connected to reservoirs, 
as there will be no possible confusion below. 

3.2. Nonlocality leading to positive WL correction 

An interesting consequence of the nature of the 
weighting factors was pointed out in Ref. [60] : since 
the weight dGaff^^/dk may change in sign for cer¬ 
tain wires, the WL correction to some transmission 
coefficient may become positive. Such an example 
is shown on Fig. [^ which has 

ir„ = l(-i + iv,^), (27) 

valid for Na long wires (^ la, h)- This WL correc¬ 
tion may become positive for sufficiently large Na, 
as a striking illustration of the nonlocality of quan¬ 
tum transport. 
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Figure 4: For a sufficient large number Na of long wires, 
the WL correction to the transmission AT 12 may he become 
positive m- 


3 . 3 . Effect of nonlocality on the four-terminal re¬ 
sistances 


We have demonstrated that a relation similar to 
(261 holds for the 4TRs, Eq. ([^ (cf. Appendix B). 
We illustrate this formula by considering the resis¬ 
tances TZi 2 ,i 2 and 7?,i2,34 characterizing the network 
of Fig. 


A remarkable consequence of quantum nonlocal¬ 
ity is the possibility of large WL correction to the 
conductance, induced by the presence of long ID 
contacts m- To illustrate this idea, we analyze 
the WL correction to the two resistances TZi 2 ,i 2 and 
"7^12,34 for the network of Fig. 

In the following we will express the resistances in 
units of hf (2^6^). 


3.3.1. Reminder : two-terminal configuration 
At this stage it is useful to recall the well- 
known result for WL correction to the conduc¬ 
tance (or resistance) in a two-terminal measure¬ 
ment (a wire of length L between two large 
contacts). In this case the dimensionless re¬ 
sistance is simply 7?.ciass = L/^ioc and the 
Cooperon satisfying Dirichlet boundary condi¬ 
tion reads Pc{x,x') = L^smh.{x^/L^)smh{{L — 
X:>)/Lip)/sinh{L/Lip), where a:< = min(a:,a::') and 
a;> = max(a:,a;'). The application of 0 gives the 
result of Ref. [B] : 


An 

^^class 


Lip 

~i7 


coth{L/L^) 


L(p 

X 



(28) 

where the expansion corresponds to <C L. The 
dominant term corresponds to the value of the 
Cooperon in bulk (i.e. inside a wire, at distance 


3> Lip from the boundaries) ; the correction is ex¬ 
plained by the depletion of the Cooperon which 
vanishes at the boundaries. In the coherent limit 
{Lp —>• 00 ), one get the universal WL correction 
An/nh^^ ~ 1 / 3 . Note that (281 may be obtained 
more directly by using the spectral determinant ap¬ 
proach [TIEB]. 


3.3.2. Two-terminal resistance 7 ^ 12.12 for multicon- 
nected wire 

We consider the device of Fig. and consider first 
the two-terminal resistance. Due to nonlocality, al¬ 
though the two arms d and / play no role for the 
classical transport, nf^f 2 = (^a + 4 + 4)/6oc, they 
influence the quantum contributions to the trans¬ 
port coefficients. Eq. 0 implies that integration of 
the Cooperon runs over the three wires (a), {b) and 
(c) : 


A7^l2,12 f . f 


/ \dxPc{x,x). 


. 

where L — la + lb P Ic- The explicit expression of 
the integral dx Pc(x, x) over a wire in terms of 
the network properties can be found in Ref. |60] . 
Below, we only analyze limiting behaviours. 


Weakly coherent limit Lp <C la, h, Ic, Id, If —• We 
get in this case 


A7?.i2,i2 ^ Lp 2 / Lp\ 


up to exponentially small corrections. If we com¬ 
pare this expression with ( [2^ , we remark that the 
presence of the two long wire (d) and (/) only af¬ 
fects the subleading term of the WL correction. The 
difference correction is explained by the depletion 
of the Cooperon at the vertices, which behaves as 
Pc{x,x) ~ Lpl2 — {Lpld)er'^^l^f at distance x of 
the vertex. The correction to the bulk result is —4 x 
(2X) dx (L^/6)e-2-/^- = -(2/3)(L^/L)2. 


Fully coherent limit (Lp -^ 00 ) —. We only dis¬ 
cuss two limiting cases as the general expression 
is rather cumbersome. When the two wires (d) 
and (/) are very long. Id, If L, h, h, they play 
no role and we recover the universal result for the 
two-terminal wire [limit Lp —>■ 00 of Eq. (28)] 

A7^12.12/(7^^'2P2)' - 1/3. 

In the opposite limit Id, If PI la, h, Ic, the phase 
coherence is broken at the level of the vertices due 
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to the vicinity of the large contacts. We get the 
result A 7 ^l 2 ,l 2 /( 7 ^^'J:f 2)2 cs ( 1 / 3 )(Z 2 +/2 + + 

lb + IcY, which corresponds to the addition of re¬ 
sistances A 7 ?,i 2 ,i 2 — + ^Rb + ^Rc for three 

independent wires with Ai?^ = (l/3)(Za/6oc)^, etc. 

Varying the length of the wires (d) and (/) hence 
allows to cross over between the fully quantum 
regime where the wires (a) + (b) + (c) should be con¬ 
sidered as a whole and the regime where the three 
wires are independent and their resistances may be 
added according to the classical Kirchhoff law. This 
is an illustration of the idea introduced by Biittiker 
to describe dephasing in a fully coherent system by 
introducing fictitious voltage probes m- 

3.3.3. Four-terminal resistance TZi 2 ^zi : nonlocality 
leading to large WL correction 

Classical resistance reads = 4/6oc which 

immediately shows that the WL correction is given 
by an integral of the Cooperon over the wire (b) 
only : 


A7^12,34 




/ da; Pc( 2 ;, x) . 

/ wire (6) 


(31) 



Figure 5: A four-terminal device. When the central wire ( 6 ) 
is shorter than the phase coherence length and the connect¬ 
ing wires, the WL correction A 7 ^i 2 , 34 /( 7 ^ 12 ^ 3 ^ 4 )^ is large. 


Incoherent connecting wires <C ^c: ^d^ — 

. Using again the expression of the integral of the 
Cooperon integrated in a wire m. we obtain the 
explicit expression 


A'Pi 2,34 -1 + Vlk coth y/^lb 




111 

-1-1- yFflb COthy^4 


(32) 


4 coth ^Ib 5 


sinh"' yFflb 


( coth y/^lb + 2 ) ( coth y/jlb 




sinh ^Ib 


where 7 = l/L^- The first term is the result for 


an isolated wire of length 4, Eq. (28). The second 


term originates from the non vanishing value of the 
Cooperon at the two vertices, i.e. we can interpret 
this term as coming from the modification of the 
boundary conditions for the wire ( 6 ) induced by the 
presence of the connecting wires (a), (c), (iP and 
(/). After a little bit of algebra, we obtain :|j 


^^ 12,34 ^ Ly, 5 coth.{lb/Lip) 4 : — 3 Lip/lb 
~ 4coth(4/Tv^) + 5 


■ (33) 


• For a long wire L^p <C 4 we get the small correc¬ 
tion 


An 


12,34 




4 


4 


(34) 


up to exponentially small corrections. The dom¬ 
inant term L^/lb <C 1 , coincides with the two- 
terminal measurement for the wire of length lb., 


Eq. (28). As for APi 2 ,i 2 the presence of the con¬ 


necting wires manifests itself only through the co¬ 


efficient of the subleading term {Lp,/lb)‘^, cf. (28) or 
(30). The subleading correction is half of the one 
obtained for APi 2 ,i 2 , Eq. (1^, as only the deple¬ 


tion of the Cooperon at the two boundaries of the 
wire ( 6 ) contributes. 

• For a short (coherent) wire 4 L,p, quite re¬ 
markably, we obtain a large WL correction to the 
four-terminal conductance 


An 


12,34 


IT. 


2 lb 


»1 . 


(35) 


In this case the presence of the connecting wires 
strongly affects the Cooperon inside the wire, which 
is the reason for the large WL correction. In a 
wire connected to two large reservoirs, the size of 
the electronic trajectories contributing to the WL 
is bounded by the length of the wire, which leads to 
a saturation of Ag as L^, —>■ 00 . On the contrary, in 
the four-terminal configuration, the electronic tra¬ 
jectories can explore the connecting wires (a), (c), 
(d) and (/) on large scales compared to lb, which is 
the physical origin for the large WL. 


Eq. (35) characterizes the WL correction to the 


conductance, in a four probe configuration. Al¬ 
though the correction to the conductance may be 


^ This result was obtained by Santhanam m, although 
this paper does not provide a detailed discussion on how the 
Cooperon must be integrated in a complex geometry. 
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large, we remark that the relative correction, 
A7^12,34 


'TP class 
'^ 12,34 


26 o 


< 1 


(36) 


is always small, as the validity of the perturbative 
treatment requires L^p ^ .fioc. 

In order to better understand the difference be¬ 


tween the two results (34) and (351, it is instructive 


to consider the multiterminal network of Fig. 
The WL correction to 7?.i2,34 is given by @ as 
well. For the calculation, the key point is that, 
at a vertex x from which issue rua, long wires 
(longer than L^), the value of the Cooperon is 
Pc{x,x) ~ L^jrax (cf. Ref. [35] and Appendix C). 
In the weakly coherent wire limit (L^p ^ 4)j the 
Cooperon inside the wire (6) is Pc{x^ x) ~ Ti^/2, ex¬ 
cept at a distance ^ from the vertices. Integra¬ 
tion of the Cooperon leads to A7?.i2,34/(7?.fj(34)^ — 
Lp/lh <C 1, similar to Eq. (341. In the coherent 



2K 


Figure 6: A wire of length connected to 2Na reservoirs by 
long wires. WL correction A7^i2 ,34/(7^fi(34)" is controlled 
by electronic trajectories starting from wire (b) which may 
explore the wires over long distances compared to If,. 


limit [Lp, <C 4)) with connecting wires still longer 
than Lip, the Cooperon is almost uniform inside the 
wire, with a value Pc{x, x) ~ L^/(2A4), where 2Na 
is the effective coordination number. We get in this 
case the large WL correction A7?.i2,34/(7?.5'j(|4)^ ~ 
Lp/{Nalb), which reduces to (35) for Na = 2. This 
argument provides the interpretation of the factor 
1/2 in Eq. (1^. 


Fully coherent limit —>■ oo —. Coming back to 
the simple network of Fig. it is also interesting 
to consider the fully coherent limit. Using again 
the expression of the integral dx Pc(x, x) given 
in |60|, some algebra leads to 

^'^12,34 _ 1 I 2 la\\d + 4||/ + 34||(i ^c||//4 , , 

(^i^"3 + 3 la\\d + lc\\f+h ’ 


where l/4||d = 1/4 + 1/4, etc. For short con¬ 
necting wire la, Ic, Id, If h, we recover the well 
known universal result ATli 2 ^ 34 /~ 1/3 
corresponding to a coherent wire between two large 
contacts, as expected. In the other limit, 4^4 = 
Ic = Id = If, we obtain a large correction 


A7?.12,34 ^14 

"^2 4 


(38) 


similar to (351 in which Lp —> la (i.e. the cutoff lim¬ 
iting the trajectories exploring the connecting wires 
is not the phase coherence length but the distance 
4 to the reservoirs). 


The possibility for large WL correction was 
pointed out in Ref. m in the rather academic sit¬ 
uation of an isolated wire. A more precise dis¬ 
cussion was provided by Santhanam m for the 
case we have considered here. As we already men¬ 
tioned, this has the same origin as the large volt¬ 
age fluctuations due to long coherent excursions of 
charge carriers in the voltage probes emphasized by 
Biittiker [16] (and also in Ref. [I3]). Although the 
observation of large WL corrections has not been 
reported so far, to the best of our knowledge, large 
resistance fluctuations have been observed in sev¬ 
eral experiments [T3l [55] . with the same physical 
origin. 


4. Fluctuations and correlations 

Mesoscopic (interference) phenomenon are more 
pronounced when the system size is reduced down 
to a size comparable to the phase coherence length 
Lp : the quantum contribution to the dimension¬ 
less conductance Ag = g — gdass of a fully coherent 
conductor presents fluctuations Sg = g — (g) ^ 1 
of the same order than the average (Ag) ^ 1 
(the WL). For this reason, the characterization 
of conductance fluctuations/correlations has at¬ 
tracted considerable attention both and experimen¬ 
tally [sansi EH EH] and theoretically p na El] 


It is also well-known that the WL correction to the 
conductance is large in the two-dimensional situation, in a 
plane [2] leading to Ag ~ —(I/tt) ln{L/£e) where L is the 
size of the plane. However, the relative correction is small, 
Agjg ~ —2ln{L/ie)/('^hFie) "C 1, where kp is the Fermi 
wave vector. The same logarithmic behaviour is obtained 
for large planar networks (square grids, honeycomb lattice, 
etc.) leading to Ag ~ —(I/tt) ln(L/a), where a is the length 
of each wire [58]. 
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(for reviews, see [Ml [71111]). The experiments are 
usually performed in the four-terminal configura¬ 
tion, which has brought the question of the role of 
the voltage probes [I3IS5|. Devices similar to the 
wire of Fig. were studied in these experiments. 
The nonlocal nature of quantum transport is par¬ 
ticularly striking by considering the symmetric and 
antisymmetric parts of the resistance 

= 2 (’^12.34 + 72 . 34 . 12 ) (39) 

= 2 (72i2,34 — 72 . 34 , 12 ) , (40) 


which were shown to present different behaviours 
as a function of the ratio Ib/L^ (symmetrisation is 
done with respect to exchange of current and volt¬ 
age probes or, thanks to Eq. (101, to magnetic field 
reversal). Whereas TZa is a relatively flat function 
of lb/Lip, the symmetric resistance presents a clear 
crossover at ^ lb, see Fig. (note that in a 
weakly disordered metal with a small enough mag¬ 
netic field, we can ignore the classical magnetore¬ 
sistance caused by the Lorentz force. As a conse¬ 
quence 72 .( 4 ®'®® = 0)- 

The study of nonlocality of voltage fluctuations 
and/or transmission probabilities in multiterminal 
devices was considered theoretically by Maekawa 
et al. |3H] and Biittiker [T^ (for a three terminal 
device) by different approaches. The importance 
of long range potential correlations was later em¬ 
phasized by Kane, Serota and Lee |43|, which has 
led to reconsider the study of voltage fluctuations 
in an illuminating paper of Kane, Lee and DiVin- 
cenzo HI], and also in Refs. |SS||37| (note also the 
numerical study El). Finally we point out few 
theoretical works on the ring configuration [HI [28] 
relevant for the experiments aforementioned, and 
specifically studied in Ref. [15] . 

In the following we derive formulae analogous to 
ffi for the correlations of transmissions and the 
correlations of four-terminal resistances 


. (41) 

We will apply our results to the analysis of the resis¬ 
tances 72-5 and TZa characterizing the four-terminal 
wire of Fig. We will show that our formalism 
allows to recover the results of Refs. m IM1I37] 
straightforwardly. 


/.I. Conductance correlations 

Expression of the conductivity correlations in 
simple geometries can be found at several places [HI 


[7]. In networks, the correlations of transmission co¬ 
efficients {5Taii{B)5Tp_v{B')) (i.e. conductance ma¬ 
trix elements) are given by four contributions [64j : 


^ f dujSxiuj) I dxdx' 

J ./Network 

dPd{a,x) dPdiii,x) 

dx dx [x,x )r_^[x,x) 

dPd{x',P) dPd{x',v) 


dx' 


dx' 


(42) 


=-4 f da;(5T(w) f dxdx' 

J ./Network 


dPd{a,x) dPdix,n) 

dx dx ^ Kx,x ,x) 

dPd{p,x') dPd{x',/3) 


dx' 


dx' 


(43) 


= ^ [ da;(5T(w) f dxdx' 

J ./Network 


Network 

PL'^\x,x')P/;/^\x' ,x) 


dPd{a,x)dPdix,(3) 

Sx - 

dPd[p,x') dPd{x',v) 

dx' dx' ^ ’ 

= same as (45) 

with ^ P(/) , 


where we used the same notation as before, 
Pd{<A,x), in order to designate the Diffuson mea¬ 
sured at a distance of the vertex a. The function 
is a normalised function]^ of width Aw ~ T 
with 5 t( 0) = l/(6r). Several remarks : 

• As for the WL, note that these expressions are 
of great generality, and not only valid for net¬ 
works of quasi-ID wires ; they only assume 
contacts of quasi-lD nature. For a more gen¬ 
eral situation, one has to replace the deriva¬ 
tives as d:^Pd{a,x)d:^Pd{p,x) -> VPd{a,r) ■ 
VPd(M,r1. 


• The Diffuson Pd{x,x') connecting the contacts 
to the bulk, and providing the weights to at¬ 
tribute to each wire, obeys the classical diffu¬ 
sion equation (13). The Cooperon and Diffu¬ 
son p!;j^’^\x,x'), which describe phase coher¬ 
ent properties, are solutions of the diffusion 


® Its precise form for is ( 5 t(w) = F(w/2T)/2r with 
F{x) = (xcotha; — 1)/sinh^ x. 
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Appendix B.4) : 


equation 


L% D ^ 


Pi‘^^‘'Hx,x') = S{x - a;') 


(46) 

where the covariant derivative = dx — 
2ieA±(a;) involves the vector potential A± = 
[A ± A']/2 for Diffuson (A_) and Cooperon 


(^+)- 


• The penetration of the magnetic field in the 
wire is taken into account through the substi¬ 
tution l/L^ 1 / 7^2 _|_ as for the 

WL (§1^. 


• These expressions are based on the current 
conserving expressions for the conductivity 
correlations, given by the procedure of Kane, 
Serota and Lee (i.e. we combine results of Refs. 
[S] and [13]). 


• Transmission correlations are related to con¬ 
ductivity tensor correlations, aab = e^voDab- 
The two first contributions ( 42|43 ), which cor¬ 
relate the indices of the two tranmissions are 
interpreted as diffusion constant correlations 
{e‘^i'o)‘^{SDabSDcd), while the two last contri¬ 
butions ( 44|45 1, which do not correlate indices, 
are related to density of states fluctuations 
{e^DY{5v'^) |3|. 


• The products of Diffusons may be related to 
derivative of classical transport coefficients, as 
it was done for the WL. For example 

1 dPd{a,x) dPd{ii,x) 1 

il dx dx (loc dli 

when X belongs to the wire i, etc. 


4-2. Four-terminal resistance correlations 


We now simplify the above expressions by ne¬ 
glecting the effect of thermal smearing, 6 t{oj) —>■ 
S{u;). Thermal effect will be described later in 
§ |4.3.3[ We deduce the correlations for the 4TRs 


(cf. 


(48) 


Sloe 

class class 




dh 


dir 


Sloe 

n' 7 ^ class class 

^''^0-13,r-'v' 


dx / dx' Pd{x,x') 
/p) J{r) 

(49) 


dh 


J Class p p 

'' / da: / da:'Pc(a:,a:')^ 
Oil' J{t) J{i>) 


{P-ap^r-vP-a'p'= JT~ 


(50) 


loc 


/ii'Tp class class 


iSS p p 

dx dx' Pd{x,x'Y 
h' J(i) J(i') 


dh dh 

p'^ — same as (***)^ ^ 

with Pd^ Pc, (51) 


where ^loc = OidNcte is the localisation length for 
the infinitely long wire. As there will be no pos¬ 
sible confusion, we now adopt the simpler nota¬ 
tion Pd^c = Pi'Z'Q . Note that, as for the transmis¬ 
sion correlations, only the contributions (• • and 
(•••)^^^ correlate the indices in a non-trivial way. 
As for the WL, the contributions of each wires are 
weighted by classical quantities. 


4-3. Four-terminal resistances in a multiterminal 
wire 


We apply our formalism to the analysis of the 
4TR correlations in a wire connected to several volt¬ 
age probes, like the one represented in Fig. Let 
us first recall the expressions of the classical resis¬ 
tances, that will be needed to compute the weights 


Eq. ( |48|49|50|51D : 


= Ib/^ioc (52) 

P12P2 = (la +lb + /c)/6oc (53) 

Pf474 = (Id + lb + h)/Zloc, (54) 


where ^loc is the localisation length in the infinitely 
long wire. As they were considered in Ref. [T3], we 
will study the symmetric and antisymmetic resis¬ 
tances TZs,a, which exhibit remarkable behaviours. 
The relations between correlators are : 

(SPIa) = I ((<57^?2,34) ± (<57^l2,34^7^34.12)) ■ 

(55) 
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We will first consider the limit of strong magnetic 
field (Lg ^ Lip or 4), when the Cooperon contribu¬ 
tions (49|51 1 are suppressed. The effect of a small 
field will be discussed in § |4.3.4| 

4-3.1. Weakly coherent regime h 

Fluctuations (i 57 ^i 2 34 ) —• The contribution (48) 
explicitly : 


IS 


A /i'Tp class class 

2 \(1) _ 4 ^'^34,34 


Cl — 


X j dx I dx' [Pd{x,x')]^ . (56) 

Fi) hi) 

The wire weights are all equal to 1/Cioc 
ply X € (a), ( 6 ), (c) and x' S (d), ( 6 ), {/), thus 
this contribution can be splitted in nine terms Ci = 
Ci.i +C 1.9 of four different types ( 57|59|62j^ ). 
The first term involves a double integral in the cen¬ 
tral wire 


leading to|^ 




( 61 ) 


Two terms involve integration in neighbouring long 
connecting wires, like 


Ci,( 


(62) 


When Lp lb, we have Ci _6 = Ci^y ~ Ci_ 2 - The 
two last contributions are of the kind 


Cl,I 


b I 

d /lx’ 


(63) 


and are exponentially suppressed 


''la 


(57) 


Ci,i = ^ / dx / dx' [Pdix,x')]‘^ 

Moc {b) J (6) 

^ ^ 4 [ d(x - x') [Pdix, x')f ~ 

Moc 


(58) 


L^h 

^loc 


where we have used the expression Pd{x,x') ~ 
(L,^/ 2 )exp [ — |x — x'\/Lip\ valid in the bulk (i.e. 
the expression obtained in an infinitely long wire), 
as the presence of the wires (a), (c), (d) and (/) 
only affects the diffuson at distance < from the 
two vertices. 

Then four terms, Ci _2 to Ci^s, involve integration 
in neighbouring wires with one coordinate in the 
wire ( 6 ) and the other in a connecting wire, like 


C 1.8 = Cl ,9 ~ 


(64) 


The contribution (< 57 ^ 12 , 34 )^^^ is simpler to dis¬ 
cuss as it involves the weight 


dh dU 


(65) 


leading to the constraint x, x' S (5). Thus we im¬ 
mediately get 


(<57^L.34)^'^ = ^Ci.i . 


( 66 ) 


Summing the ten terms, we get the fluctuations 
((577.22 34 ) ” (^/^) Ci,i 4“ Cl, 2 -f *'' Cl, 9 . Since the 
first term dominates, we recover the well known 
expression [5] for the conductance fluctuations 


■' 1,2 


(59) 


In the limit L^ ^ lb, integrals are dominated by 
X and x' close to the vertex and we can use the 
expression of the Diffuson for the star graph with 
coordination number 3 given in [Appendix C| 


Cl ,2 - jr- / 

^loc 0 


dx / dx' f 


(60) 


(( 577 . 12 , 34 ) d f Lp 


2 V k 


X « 1 • 


(67) 


Correlations (d 77 i 2 , 34 ^ 7734 , 12 ) —. The weights of 
the contribution (• • are 


dTZfYii dTZfYi2 

dk dL 


( 68 ) 


® The equality Ci ,2 = ■ ■ 
wires much longer than L^. 


= Cl ,5 holds for connecting 
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which requires once again x, x' S (6). Therefore 

(<572.12,34^72.34,12)^^^ =Ci,i- 

The contribution (•••)^^^ takes the same form, 
hence (572.12,34(572,34_i2)*'^^ = (l/2)Ci,i. 

The correlations are given by summing these two 
terms ( 572 . 12 , 34572 . 34 , 12 ) = (3/2)Ci,i, which coin¬ 
cides exactly with the dominant term of the fluc¬ 
tuations. 


Conclusion —. If we now gather all these results 
we obtain 


(<57^1) 

{sn\) 


3Ci i -I-Ci 2 + ■ ■ ■ + C 


1,9 


Cl,2 + 


2 

+ Cl ,9 


3T^4 

2^ 

4 


^loc 


(69) 

(70) 


6TZs grows with 4 whereas STZa is independent on 
the distance between the two voltage probes, in 
agreement with the experiment (see figure]^. 


4-3.2. Coherent limit l^ 

We now discuss the case where the central 
coherent, with long connecting wires 


wire IS 


(4, 4j h, If iS> L^). We can write that the Dif- 
fuson is almost uniform inside the wire (6), equal 
to the Diffuson at a vertex of coordination number 
4. Thus Pd{x,x') ~ T(^/4 for x, x' € (5). It decays 
exponentially over the distance in the connect¬ 
ing wires : Pd{x,x') ~ [L^/^) exp[-(x -I- x')/L^p\ 


for X S (a) and x' S (c), etc. (cf. Appendix C). As 
a result we see that the fluctuations are dominated 
by four terms, Ci,6, • • • ,Ci ,9 corresponding to cases 
like (621 or ( [^ , when x and x' are both integrated 
over the long distance in the connecting wires : 
( 572 .^ 2 , 34 ) — (1/4) (L,^/4oc)^, i-e. the four-terminal 
conductance fluctuations are large 


( 57212 , 34 ) 


=.4171 »i 


(U) 


The correlations are much smaller as they involve 
integrations of x and x' in the central wire only 


( 572 . 12 , 34572 . 34 , 12 ) ~ (3/8) T^^b/4oci i-6- 


( 572 . 12 , 34572 . 34 , 12 ) ^ 3 
- 8 

Finally we get 


> 1 


(57^|) C4 (57^i) C4 ^ (^1^ 


(72) 


(73) 


Fluctuations are independent on the distance lb be¬ 
tween the voltage probes. 


4 . 3 . 3 . Thermal smearing 

In the two previous paragraphs, we have repro¬ 
duced the main conclusions of Hershfield [37] by 
simpler arguments based on the analysis of the wire 
weights. Although the qualitative change of be¬ 
haviour between STZs and 6 TZa at lb agrees 

with the experiment of Benoit et al. [13] (Fig. in, it 
was noticed that thermal smearing, not accounted 
for by Hershfield, is important in the experiment. 
This corresponds to the case where the thermal 
length Lt = D/T is smaller than L^. 

There is also a fundamental reason to consider 
this regime : in the low temperature regime (T < 
1 K) and in the absence of external degrees of free¬ 
dom like magnetic impurities, the decoherence is 
dominated by electronic interactions. As a conse¬ 
quence one has Lt ^ L^ (see [T] [SIl HI [SS] |S5] and 
references therein). 

/o'! 

In this regime the contribution {■ ■ -y is negligi¬ 
ble compared to (• • •)^^^ 


Weakly coherent limit Lt ^ L^ lb —. We 
should repeat the analysis of the previous subsec¬ 
tion by adding an integration over the frequency 
with the thermal function. The first contribution 
to the resistance fluctuations reads 

Ci,i = ^ f dujST{uj) (74) 

X f dx f dx' P^‘^\x,x')p!:^^{x,x') . 

J(b) J{b) 

At finite temperature, the Diffuson involves expo¬ 
nentials like exp[— 77 a;] with 7 = 1/T^ — iuj/D 
where uj/D < 1/L|,. When min(LT,T<^) ^ lb, the 
Diffuson decays rapidly inside the wire. Ci,i is dom¬ 
inated by the integral in the bulk of the wire, which 
leads to a similar calculation as for an infinite wire. 
When Lt «C L^ (the case L^ Lt lb is similar 
to Lt = 00 , and was discussed above), the ther¬ 
mal function can be considered as a broad function 
and simply replaced by 5t(w) —>■ I/( 6 T). Using the 
expression of the Diffuson in bulk we get 

Cu 4 , (75, 

Sloe 


This term dominates the fluctuations, therefore we 
recover the known expression [5] of the conductance 
fluctuations 


(572.^2,34) ^ 7 ’’ L'^Lp 


(76) 
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Next contributions are the six terms of type (59) 


and (62). We can again argue that integration is 


dominated by x and x' close to the vertex and use 
the approximate form given in [Appendix C| : 


r 4 /• du; 

“ Cfoc ./ 


(77) 


X / dx f dx'P^'^\x,x')P^^{x,x') . 

J(a) J(b) 


where Pi, (a;,x') ~ l/{3y/j)e ) with 7 = 

1/L^ — iujjD. Some algebra gives 


Cl 2 — ■ ■ ■ — C-] 


2 L'lL 


1.7 




27 ^,1 


(78) 


Correlations are obtained by similar arguments 
(<572-12,34^7?,34, 12 ) = C 17 , which coincides with the 
dominant term of the fluctuations. 

Going back to the quantity of interest, we should 
replace 
whereas 


= (79) 

Sloe 

(80) 

^ Sloe 


^by { 6 nl) = Ci,i - 
(|70|) still holds, thus 


(^ 1.2 


+ Cl, 9)72 


We reproduce similar conclusions as for Lt = 00 , 
i.e. fluctuations STZs growing with the distance lb 
between the voltage probes, while the fluctuations 
STZa are independent on lb- 


Coherent limit lb At ^ L^p —. As in the case 
lb Lip <ti. Lt analyzed previously, the fluctu¬ 
ations are dominated by four terms where x and 
x' are integrated over two long connecting wires. 
It is simply given by multiplying (78) by a factor 
(3/4)^, which accounts for the fact that the Diffu- 
son feels an effective coordination number 4 instead 
of 3, hence ~ ~ C^g ~ (1/24)L|L2 

Thus 


in order to cut off the contribution of large fre¬ 
quencies. When Lt <C L^p, some algebra gives 
Ci,i (l/6)(Af./^/5f„Jln(L,^/LT), leading to 


(< 572 . 12 , 34 ^ 72 . 34 , 12 ) 

class '77 class \2 
V^12,34 '^34,12.' 


1 

6 



2 

ln(A,^/AT) ■ 


(82) 


Quite surprisingly, we obtain a logarithmic depen¬ 
dence in Lip reminiscent of the 2D situation, al¬ 
though the system is ID. 

The symmetric and antisymmetric resistance 
fluctuations are 

(57^|)=:<(^7^i):x 1^1^. (83) 

Sloe 


4.3.4- Magnetic field dependence 

Finally, we discuss the magnetic field depen¬ 
dence. The above results are valid when the mag¬ 
netic field B = B' is larger than the correlation field, 
i.e. in a diffusive wire, when B ^ Be ^ (po/{Lipw) 
where w is the width of the wires and fio = hje the 
flux quantum. At small magnetic field, Cooperon 
contributions ( 49|51 ) must be taken into account 
as well. This can be done easily by noticing that, 
due to the wire weights, the role of fluctuations and 
correlations are exchanged for the two first contri¬ 
butions : 

(<57^?2,34)^"^ = (^7^l2.34<57^34,12)^'^|p^^p^ 
((57^l2,34<57^34,12)^"^ = (^7^?2.34)^'^|p,^p^ 
while 

(<57^?2,34)^"^ = (<57^?2,34)^'^ 


(57^l2.34<57^34,12)^"^ = (<57^ 


The substitution Pd 
doing 1/Ll l/T'v 


Pd^Po 

1^72..s4.12)^^^ 


12,340/V34^12; 


Pd—^Pc 


Pr is simply achieved by 

2 cf. § 


1/7'B) 


3.1 


(the anal¬ 


ysis of the correlations for different heids requires 
similar substitution in both the Diffuson and the 
Cooperon, cf. § 4.1). We deduce straightfor¬ 


wardly the crossover functions. We only consider 
the weakly coherent regime. Lip <C h- 


(< 572 ^ 2 , 34 ) ^ 1 [ ^ d 


(81) 


The analysis of the correlations (< 572 - 12 , 34 < 572 . 3472 ) = 
Ci.i is as follows : we write Pi‘^\x,x') ~ 

1/(4\/7) when x, x' G ( 6 ), therefore Ci,i ~ 
(l/4)(^b/fD/dw<5T(w)/l7l where 7 = l/A^ _ 
iuj/D. The presence of the function 6t{uj) is needed 


Regime Lip h, Lt —• Using the previous calcu¬ 
lations we obtain 

(«R|)-57^(i + [i + (e/s=)T"'b (84) 

^loc 

(85) 
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where BjBc = L^jL^ = {2TTl'/3)BwLipj(j)Q. 
Whereas the symmetric resistance fluctuations are 
doubled at small field, the antisymmetric resistance 
vanishes as STZa{B) ^ {L^/i\ocf {BL^w/(f)o) for 
B^O. 


Regime Lt ^ Lip R —. In this case we obtain 


{STZl)e.^^^[l+[l + {B/B^r]-^^^) 


(snl) 


3 efoc 


1 L^^Ll 


B^ 


12 ^ 4 ^ ^ 2 +^ 2 - 


( 86 ) 

(87) 


The antisymmetric resistance now vanishes at small 
field as 5TZa{B) ^ (LTLip/^f^^) {BL^w/ffio) for B 

0 . 


4-3.5. Experiments 

We now rediscuss the experiments at the light 
of our results. In Fig. [Tj we have reproduced the 
experimental data obtained by Skocpol et al. [55] 
for a Silicon inversion-layer narrow wire. We now 
denote by L the distance between the two voltage 
probes (denoted R above). The experimental result 
exhibits the behaviours obtained above : a growth 
of the fluctuations with the length, STli 2 34 ~ 
rTr'/ViL in the incoherent regime < L, 
Eq. ( [^ , and a saturation i 577 .i 2,34 ~ (Lip/^ioc)^ in 
the coherent regime L^p > L, Eq. Note that 

accounting for thermal broadening does not change 


this conclusion, cf. Eqs. (76) and 




300 



Figure 7: Mesoscopic (sample to sample) resistance fluctua- 
tions for multiconnected silicon inversion-layer narrow wire 
at T 400 mK. The two samples have lengths L ~ 0.15 pm 
r^i 2,34 ~ 450 and 5 pm (^7^i2,34 ~ 12 kQ/ Data from 
Ref. ]55f. 

Another remarkable experiment is the one of 
Benoit et al. ca, who analyzed voltage fluctuations 
in Sb and Au narrow wires. Using the data given 
in this reference, we have plotted the fluctuations 
of the symmetric and antisymmetric resistances in 


Fig.® The most striking outcome is the compar¬ 
ison of the different behaviours for the symmetric 
and antisymmetric resistances. 



Figure 8: Mesoscopic resistance fluctuations for multicon¬ 
nected Au and Sb wires at T = 40 mK and 300 mK. Length 
L varies from 0.2 pm to 4 ^m. L is the distance between 
the voltage probes. Blue circles correspond to the symmetric 
resistance and red diamonds to the antisymmet¬ 
ric resistance The inset shows the Sb sample ; 

the fluctuations for a wire of length 
Lp. Data from Ref. m- 


Let us analyze the behaviours more into details. 
We follow the line of Ref. [U] , where thermal broad¬ 
ening was not taken into account. Voltage (sample 
to sample) fluctuations were measured respectively 
to the fluctuations for a wire of length L^p, which 
coincides with the relative resistance fluctuations 


5Tls,AlSTiip where STlp = (L,^/^ioc)^- 

In the coherent limit {Lip > L), Eq. (73) gives the 
value 6TZs,a/STZip ~ l /-\/8 ~ 0.35 consistent with 
the experimental data. 

In the incoherent limit {Lip < L), Eq. (70) gives 
STZaISTZip ~ I/VS ~ 0.577, whereas (69) leads to 
STZs/STZip ~ \/3/2 yjL/Lip. Although the linear 
behaviour with ^LjL^ agrees qualitatively with 
the experimental data, the prefactor obtained ex¬ 
perimentally, « 0.47, is significantly smaller than 

VWS - 1-22. 

The fact that the experimental values are smaller 
than the theoretical predictions could be explained 
by the effect of thermal broadening, as it brings 
a reduction factor LtP^. A more precise analy¬ 
sis would nevertheless be needed. Another aspect 
which could explain this quantitative disagreement 
concerns the determination of the phase coherence 
length, which was not obtained by a unique inde¬ 
pendent procedure in Ref. [13] . At the highest tem¬ 
peratures, the phase coherence length was obtained 
by independent WL measurements, which does not 
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account for the fact that, in heavy metals like gold, 
when spin-orbit scattering is strong and magnetic 
impurities present, the phase coherence length in¬ 
volved in the weak localisation [35] and the conduc¬ 
tance fluctuations [33] differ (see also m)- The lack 
of a formula for the WL correction in the multiter¬ 
minal wire did not allow the authors to extract L^p 
for the lowest temperatures. A more reliable pro¬ 
cedure would have been to obtain from another 
long wire L^) made under the same conditions. 


5. Conclusion 


Since the pioneering work of Markus Biittiker, 
four-terminal resistance is now recognized as a 
quantity of major importance in mesoscopic physics 
(Ref. [14] is indeed Biittiker’s most cited article, 
with now almost 2000 citations). After having 
briefly reviewed several aspects of four-terminal 
resistances, we have focused the discussion on 
the analysis of quantum transport in networks of 
weakly disordered wires, for which explicit expres¬ 
sions for the quantum corrections to the conduc¬ 
tances and the four-terminal resistances were dis¬ 
cussed. We have seen that both the weak local¬ 
isation correction and the correlations involve in¬ 
tegrals of the Cooperon (and also of the Diffuson 
for the correlations) inside the wires, whose con¬ 
tributions must be weighted by derivatives of the 
classical coefficients : dG'^'^^/dk for the conduc¬ 
tances and dli for the four-terminal re¬ 

sistances. Although we had discussed this earlier 
for the WL correction to the conductance [60], we 
have provided here a new interpretation of these 
coefficients in terms of ^^generalised conductances” 
Qi^a relating the current in a wire i inside the 
network to the external potential Va, precisely : 
QQdL^ss^dli oc Qi,aQi,p- 

We have illustrated the efficiency of our formal¬ 
ism by considering simple examples. Our determi¬ 
nation of the four-terminal resistance correlations 
only involves simple calculations (which greatly 
simplify the analysis of Ref. m in particular). 
Moreover we have been able to consider the effect 
of thermal broadening, which was not studied so 
far. All the main dependences are summarized in 
Tab. [T] and Tab. [2] 



L(p If) 

Lip 


fluct. ~ correl. 

fluct. » correl. 

Lp, <C Lt 

Lt Lp, 

( 

\ h i lb 

{K/kY {Lp/kf 


Table 1: Main behaviours of the four-terminal conductance 
fluctuations 34 }/("^ 12 ^ 34 )"^ and the conductance corre¬ 

lations ((^7^12, 34 ^^ 34 , 12 }/(’^i 2 ^ 34 )‘^ of the four-terminal re¬ 
sistances. Precise numerical prefactors are given in the text. 



Lp 

< k 

Lp 



» 


Lip Lt 

/ jb_ 

Uloc/ boc 

(T^/6oc)" 

(T^/6oc)" 

Lt <C Lip 

/ Lt \ ^ lb 

Uloc7 eioceioc 

itfitf 

itfitf 


Table 2: Main behaviours of the symmetric and antisym¬ 
metric part of the four-terminal resistances, ^loc = otdNc^e 
is the localisation length of the wire (we recall that the va¬ 
lidity of the diagrammatic appraoch is ^\oc ^ min(/ 5 , L(^). 
Precise numerical prefactors are given in the text. 


Appendix A. Conductance matrix, re¬ 
sistance matrix and four- 
terminal resistances 

Appendix A.l. Conductance and resistance matri¬ 
ces 

A natural way to characterize the linear response 
of a multiterminal structure is to introduce the con¬ 
ductance matrix relating the voltages at the con¬ 
tacts to the currents : la = The ma¬ 

trix elements of the conductance matrix are not in¬ 
dependent as they must satisfy two types of con¬ 
straints : (i) current is conserved, thus J2a = 0 
whatever the choice of external potentials, (ii) A 
global shift of all voltages Ra —t 14, -I- Uq does not 
change the current. As a consequence J2a = 
^a /3 = 0. In other terms, introducing the vec¬ 
tor Xq = (1,1,-- - ,1) where denotes transposi¬ 
tion, we rewrite the second condition as GXq = 0, 
i.e. G is a linear map acting in the vector space 
£± = {X G I Xq ■ X = 0}. The current conser¬ 
vation rewrites X^G = 0, meaning that G maps £± 
onto itself. In the subspace £±, we may invert the 
conductance matrix, which leads to introduce the 
resistance matrix i? relating the external currents 
to the voltages V — RI. The relation between the 
two matrices is thus RG = GR = — T]|, where 

P|| = (1/A)AoA(]’ is the projector on the vector 
Aq, where N is the number of contacts. We write 
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more conveniently : 

'y ^ — y ^ Ga^R^X — ^qA ■ (-^*1) 

7 7 

Since we will have to consider differences between 
potentials at various contacts, the relation 

y ^ (-^a7 Rfi^') G^x — y ^ Gx^ {^R-yot Ryfi) 

7 7 

= 5aX — 5px (A.2) 

will be useful. 

Appendix A.2. Four-terminal resistance 

The four-terminal resistance provides informa¬ 
tion in the situation where current only flows 
through contacts a and (3, la = —Ip = I, with 
all other currents vanishing, Eq. ([^. The voltage 
at contact /i is thus conveniently expressed in terms 
of the resistance matrix as = (i?^a — R^.p)I and 
the voltage difference = 'R-ap,^ivI thus in¬ 

volves 

— R^a Ru'd Rva ~t“ Ri>P • (-^*^) 

A useful relation is obtained by differentiating the 
relation [ig : 

^ ^A/S) ■ 

(A.4) 

Appendix B. Averaging functions of the 
conductances 

As we have discussed in the body of the text, 
the quantum contributions to the transport coeffi¬ 
cients are naturally expressed for the eonductance 
matrix, as linear response theory provides formu¬ 
lae for the conductivity or the conductance. On 
the other hand, resistances, which are in general 
complicated functions of the set of all conductance 
matrix elements, are more easy to handle in most 
situations, and usually the relevant quantities in 
most experiment. We show in this appendix how 
one can go from the quantum contributions to the 
conductance matrix elements to the equivalent con¬ 
tributions to the four-terminal resistances. 

Let us consider a general quantity, function of 
the conductance matrix elements A = f{G), where 
f{G) is a short notation for a function of all ma¬ 
trix elements /(Gn, G 12 , G 13 , • • •) ! tti® quantity A 


may be for example a four-terminal resistance. We 
now show how the quantum contributions to the 
quantity A can be related to the quantum contri¬ 
butions to the conductance. For this purpose it 
is convenient to analyze the scaling with the num¬ 
ber of channels Nc, which is a large parameter. 
Considering the set of conductance matrix elements 
{Gap}, the disorder average is given by two contri¬ 
butions : a classical term and the weak localisa¬ 
tion correction : (Gap) = -|- AGap- Other 

quantum contributions of interest are the meso¬ 
scopic (sample to sample) fluctuations, defined by 
Gap = {Gap) + dGap and characterized by the 
correlation functions {6Gap dG■ As we have re¬ 
called in the text, the three quantities present the 
Nc dependences 


Gclaas ^ 

(B.l) 

AGap = 0{N°) 

(B.2) 

{SGap SG^A = 0{N°) . 

(B.3) 


We now deduce two useful properties. 

Appendix B.l. Property 1 

Writing Gap = {Gap) + dGap, the fluctuation 
is smaller than the average, 5Gap = 0{Nc), and 
vanishes on average by definition. As a consequence 

{f{G)) = f{{G)) + 0{fxNc-A- (B.4) 

Appendix B.2. Property 2 

Using this property, we can split the average into 
a classical part and the WL correction : (G) = 
Qciass _|_ Hence the average of the physical 

quantity takes the form 

(A) = + AG) + 0{f X A-2) (B.5) 

= /(G"") + ^ AG„„ + 0(/ X IV-^). 

/.A,I/- 

We identity the second term as the weak localisa¬ 
tion correction to A : 

n ^class 


^ For a simple geometry, the relevant large parameter is 
rather the dimensionless conductance g (which is simply g = 
^loc/-^ for ^ wire of length L, with ^loc oc Ncie)^ however in a 
complex network of metallic wires, involving several lengths, 
there is no unique such parameter. 


17 






Appendix B.3. Property 3 

If we consider two quantities A and B functions 
of the transmissions, following the same lines, we 
deduce the expression for the correlation function 


{5ASB) 


E 


^^class ^^class 

dG^dG^ 


{SG^pSGf,,) . 


(B.7) 


Appendix BA- Applications 


• WL correction to the four-terminal resis¬ 
tances : It is now straightforward to get the 
weak localisation correction to the four-terminal 
resistance. We deduce from the property (B.4|P 
- (Rm)) (G^a) = - Spx + 0{N^. 

Then expanding the con ducta nces as (G^a) = 
G^® _l_ AG.yA) and using (A.4) and the property 


(B. 6 ), one get 


ATZa 




class pclass 

/j,7 “^1^7 


(B. 8 ) 

\ A / Dclass DClaSSN 

)I\Lf^x(nxa )• 


7A 


Using ( [M| and once again the equation (A.4) we 
finally get (|^. Note that this result could have 
been guessed from the heuristic argument presented 
in the introduction of Ref. EO]. 

• Correlations of four-terminal resistances : 

The expressions for the correlation functions 
( |48|49|50|51 ) are demonstrated by making use of 
the relation ( |AA| ) with the property ( |B.7[ ) (with 
the remark closing paragraph |4.1[). 


Appendix C. Solution of the diffusion equa¬ 
tion in the star graph 

The solution of the diffusion equation [7 — 
dl]Pc{x,x') = 6{x — x') in the star graph with 
nria infinite wires is useful (Fig. |C.9[ ). Details are 
given in Appendix D of Ref. |58j . The value of 
the Cooperon at the vertex a is inversely propor¬ 
tional to the coordination number of the vertex : 


® Whereas the classical resistances and conductances sat¬ 
isfy the relation = &a\ — , the average 

resistance and conductance do not, up to correction of order 
due to the correlations between transmissions. How¬ 
ever the relation is obviously satisfied for a given realization 
of disorder : Y.p RapGpx = &a\ - IW- 


Pc{a,a) = When the two arguments 

belong to wires x G (i) and x' G (j), we get 


Pc{x,x') = - 

mayjA 


1 

'v7- 




(C.l) 


where a;< = min(x,x'), x> = max(a;,x') and the 
distance is measured from the vertex. At large dis¬ 
tance of the vertex, we recover the result for the in¬ 
finite wire Pc(ai, = [ 1 /( 2 a/ 7 )] exp [—y7|ai —a;'|]. 



Figure C.9: A star graph with coordination number = 7. 
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